Dissipation and Vortex Creation in Bose-Einstein Condensed Gases 
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We solve the Gross-Pitaevskii equation to study en- 
ergy transfer from an oscillating 'object' to a trapped Bose- 
Einstein condensate. Two regimes are found: for object veloc- 
ities below a critical value, energy is transferred by excitation 
of phonons at the motion extrema; while above the critical 
velocity, energy transfer is via vortex formation. The sec- 
ond regime corresponds to significantly enhanced heating, in 
agreement with a recent experiment. 
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The existence of a critical velocity for dissipation is 
central to the issue of superfluidity in quantum fluids. 
The concept was first introduced by Landau in his fa- 
mous criterion 0] , where elementary excitations are pro- 
duced above a velocity vl- In liquid 4 He this process 
refers to the excitation of rotons, with Vl — 58 ms -1 . 
However, much smaller critical values are observed ex- 
perimentally, which prompted Feynman to propose that 
quantized vortices may be responsible 0] . 

Vortex nucleation in superfluid 4 He is difficult to ex- 
plain quantatively. Strong interactions within the liq- 
uid, plus thermal and quantum fluctuations, impede for- 
mulation of a satisfactory microscopic theory. In con- 
trast, Bose-Einstein condensation (BEC) in trapped al- 
kali gases HH provides a relatively simple system for 
exploring superfluidity. Weakly-interacting condensates 
can be produced with a negligibly small non-condensed 
component. This allows an accurate description by 
a nonlinear Schrddinger equation, often known as the 
Gross-Pitaevskii (GP) equation. The system also offers 
excellent control over the temperature, number of atoms 
and interaction strength, as well as allowing manipula- 
tion of the condensate using magnetic and optical forces 
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Recent experiments have produced vortices by coher- 
ent excitation Q and cooling of a rotating cloud |J . The 
existence of vortices was also inferred by Raman et al. 
, where a condensate was probed by an oscillating laser 
beam blue-detuned far from atomic resonance. The op- 
tical dipole force expels atoms from the region of highest 
intensity, resulting in a repulsive potential. Although 
vortices were not directly imaged, significant heating of 
the cloud was observed only above a critical velocity, indi- 
cating a transition to a dissipative regime. This heating 
was found to depend upon the existence of a conden- 
sate, indicating that it must be due to the production of 
elementary excitations, that subsequently populate the 



non-condensed fraction. 

Critical velocities for vortex formation in superflow 
past an obstacle have been studied numerically by solu- 
tion of the GP equation in a homogeneous condensate 
||-|l0]]. Simulations have also confirmed that vortices 
are nucleated when a laser beam is translated inside a 
trapped Bose condensed gas [ pdp^ ]. In this paper, we 
attempt to clarify the role of vortices in the MIT exper- 
iment by presenting 2D and 3D simulations of an os- 
cillating repulsive potential in a condensate. The motion 
transfers energy to the condensate, and we observe that 
the transfer rate increases significantly above the critical 
velocity for vortex formation. 

Our simulations employ the GP equation for the 
condensate wavefunction ^r^t) in a harmonic trap 
Vtrap(r') = Y^j^jJ 2 '' = x,y,z. For convenience, we 
scale in harmonic oscillator units (h.o.u.), where the units 
of length, time, and energy are (h/2muj x ) 1 ^ 2 , lu^ 1 , and 
fko x , respectively. The scaled GP equation is then 



idt® = (-V 2 + V + C|*| 2 )*, 



(1) 



where V represents a time-dependent 'object' potential 
superimposed upon a stationary trap: V = \{x 2 + ey 2 + 
r/z 2 ) + V b{f,i)- The atomic interactions are parameter- 
ized by C = {NUo/hw x ){2muj x /h)' 1 / 2 : N atoms of mass 
m interact with a s-wave scattering length a, such that 
?7o = Ai:h 2 a/m. The number of dimensions is 7. For 
most of the simulations here, 7 = 2, corresponding to 
the limit r\ — > 0. In this case, N represents the number 
of atoms per unit length along z. For the general 3D 
situation, a moving laser beam focused to a waist wq (in 
h.o.u.) at (0, y'(t),0), is simulated using 



V oh (r,t) = — exp 
a 



-2(x 2 + (y y'(t)) 2 ) 



aw. 



7,2 



(2) 



where a = 1 + (z/zq) 2 . The Rayleigh range is zq = 
ttWq/X, where A is the laser wavelength jL3|. 

Our numerical methods are discussed elsewhere . 
Briefly, initial states are found by propagating (Q) in 
imaginary time with V^bt^jO), using a spectral method. 
Then, real-time simulations are performed subject to mo- 
tion of the object potential. To recover the essential 
physics behind the MIT experiment @, we describe the 
oscillatory motion by y'(T) = a — vT (T < 1/2/) and 
y '(T) =vT-3a (1/2/ < T < 1//), where T = t — s/f 
and s is the number of completed oscillations. The veloc- 
ity between the motion extrema is constant, v = ±Aafy, 
where a is the amplitude and / is the frequency. The con- 
densate is anisotropic, with its long axis along y (e < 1). 
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As a consequence, for small a, the beam moves through 
regions of near-constant density. Initially, the object cre- 
ates a density minimum at y — a, which follows closely 
behind the moving object. For v > v c , where v c oc c s 
and c s — y/2C\ x &\ 2 is the sound velocity, the density in- 
side the beam evolves to zero. This is accompanied by 
a 7r phase slip I jUJ , at which point the density minimum 
splits into a pair of vortex lines of equal but opposite cir- 
culation |L5[ . The vortex pair separates, and the process 
begins again. 

The creation of phonons or vortices increases the en- 
ergy of the condensate, which was calculated numerically 
using the functional E = J (|V*| 2 + V\^\ 2 + §|*| 4 )<2 3 r. 
The time-independent ground state of the wavefunction 
represents the minimum of this functional. The energy is 
related to the drag force on the object .F b by dE/dt = 
F a b.v. The drag can be calculated independently over 
the whole condensate using F b — — J I^^VV^bC? 3 ^, al- 
lowing a numerical check. Superfluidity corresponds to 
the situation where E remains constant when V b is time- 
dependent; i.e. when there is no drag on the object. 

Fig. |l| shows the energy and drag as a function of time, 
as calculated for two different frequencies in 2D simu- 
lations. At low frequency, the energy transfer is rela- 
tively small and characterized by 'jumps' at the motion 
extrema, whereas at higher frequency the energy transfer 
is two orders of magnitude larger and more continuous. 
Further insight can be gained by considering the drag. At 
low /.there is little drag except at the motion extrema 
(Fig. 0(c)), while at high / appreciable drag is observed 
at all times (Fig. g(d)). 
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FIG. 1. Time-dependent 2D simulations of laser beam os- 
cillation, with grid spacing of 0.156 (512 x 128 points) and 
parameters C = 1000, e = 0.0625, a = 4, U oh = 20 and 
wo — 1.0. Condensate energy as a function of time are plot- 
ted for (a) / = 0.05 and (b) / = 0.2. The drag F oh is also 
plotted for both frequencies in (c) and (d) respectively. 

To measure the average rate of energy transfer, a linear 
regression analysis is performed on the energy-time data. 



The gradients are plotted against v in Fig. 0. It can be 
seen that the curves are characterized by two different 
regimes. Small energy transfer at low v, gives way to 
enhanced heating above the critical velocity, v c . At high 
v, the three plots follow a single linear curve. 
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FIG. 2. Mean rate of energy change as a function of ve- 
locity, for a = 3 (triangles), a = 4 (squares) and a = 5 
(bullets). Otherwise, parameters are the same as Fig. |IJ The 
dashed line shows the speed of sound in the condensate center, 
c s = ^/2jl ~ 3.55. The plot shows a sharp transition between 
phonon heating (low v) and vortex heating at v c ~ 0.4c s . 

Energy transfer below v c arises due to emission of 
sound waves at the motion extrema. This process (hence- 
forth referred to as phonon heating) is found to approx- 
imately scale with u 3 , indicating that at each extremum 
(which are reached at a rate cx v) a sound wave with en- 
ergy ~ v 2 is emitted. Note phonon emission by the object 
is not inconsistent with Landau's criterion. In particular, 
the Landau argument relies on use of Galilean invariance, 
which breaks down when the condensate density varies, 
or when the velocity changes abruptly. 

For the parameters we have explored, phonon heating 
is found to be relatively small compared to the energy 
transfer from vortex formation above v c . The heating 
rate in the latter regime is found to scale approximately 
linearly with v. This implies that the drag force is con- 
stant. Indeed, we observe that the drag saturates as v 
increases. This behavior contrasts with that of steady 
flow, where the drag oc v (where k ~ 1 at v close to 
v c , and k — > 2 for v > c s ) ||[l0|. The difference arises 
from the oscillatory motion: as the object travels back 
through its own wake, a large pressure imbalance across 
the object does not develop. 

Fig. H plots the mean energy transferred against the 
number of vortex pairs (counted in the simulated wave- 
function) . The energy transfer per vortex pair is approx- 
imately constant, leading to an estimate of the pair en- 
ergy, which is plotted as a function of nonlinearity in Fig. 
| (inset). The energy of a vortex pair in an homogeneous 
condensate with number density n, is given by 
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where £ is the healing length and d is the distance be- 
tween the vortices. Equation (^) is valid for the inhomo- 
geneous condensate when £ <C d <C R, where R is the 
radial extent of the condensate. Eq. (||) with d — 2u>q is 
plotted in Fig. || (inset), and is found to agree with the 
numerical data. Recall that the vortex pair separates 
immediately after formation, when the pair still resides 
within the density minimum created by the object. The 
pair also moves in the direction of the object motion: 
however, it is slower, and is eventually left behind. At 
this point, it has an energy approximately equal to -E pa ir, 
and the formation process is complete. The heating rate 
can be expressed as dE/dt = E pa i T f s 0, where f s is the 
shedding frequency, which is found to be proportional to 
v. This accounts for the linear dependence of the energy 
transfer rate. 
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FIG. 3. Number of vortex pairs created up to t = 10 
against rate of energy transfer. Simulation parameters are 
the same as Fig. [j], with a = 4 and C = 2000 (plotted with 
squares, fit with a solid linear regression line); C = 3000 
(circles, dotted line); C = 4000 (diamonds, dashed); and 
C — 5000 (triangles, long-dashed). The data points closely 
follow the regression lines, suggesting a constant energy for 
each vortex pair. Inset: the average pair energy against C, 
where the dashed line shows the pair energy predicted by Q . 

The subsequent vortex dynamics involve an interplay 
between velocity fields induced by other vortices, and ef- 
fects arising from the condensate inhomogeneity. In the 
absence of the object, an isolated pair follows a trajectory 
similar in character to that of a vortex ring Jl4[ , culmi- 
nating in self-annihilation. However, the object moves 
back through its wake, interacting with the original pairs 
and creating more vortices. The circulation of a pair de- 
pends upon the direction of the object motion when it 
is created. So, vortex pairs of opposite circulation are 
formed and interact when sufficiently close. This leads 
to situations where vortices annihilate or move towards 



the edge. The number of vortices remaining within the 
condensate bulk is found to reach an equilibrium value. 
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FIG. 4. Critical velocity for vortex formation v c at 
C = 2000 as a function of potential height U a b, expressed 
as a fraction of the speed of sound at the condensate center, 
c s . Inset: critical velocity plotted against C, with U b ~ 20. 
The other parameters in both plots are a = 4 and Wo = 1. 

The critical velocity for vortex formation, v c , as a func- 
tion of potential height and nonlinear coefficient is shown 
in Fig. [|. The critical velocity is not as well defined as 
in the homogeneous case [p|-^0[ for a number of reasons. 
First, a density inhomogeneity along the direction of mo- 
tion leads to a variation in c s , and therefore v c . However, 
this is less than ~ 3% in the simulations considered here. 
The oscillatory nature of the object motion is important. 
The time taken for a vortex pair to form diverges to in- 
finity as v approaches v c from above. So, the measured 
value of v c increases from its true value as a decreases. In 
addition, the object travels through its own low-density 
wake, where c s is lower. Vortices can therefore be formed 
after the first half-oscillation, when v is slightly below 
v c . Nevertheless, we can obtain a good estimate for v c 
by choosing intermediate values of the amplitude (e.g. 
a = 4) and considering only vortex formation during the 
first half-cycle. 

Fig. [I| demonstrates that v c decreases as a function 
of increasing object potential height, U a h, allowing an 
experimental diagnostic for vortex formation at varying 
beam intensities. This behavior agrees with simulations 
of ID soliton creation |l6| and vortex ring formation in 
3D p7| . We have also studied the case of E/ b < 0, which 
corresponds to a red-detuned laser. Atoms are attracted 
to the potential minimum, creating a density peak which 
moves with the beam. Vortex pairs are created from a 
density minimum which develops ahead of the beam. Fig. 
H (inset) shows v c as a function of C. The critical velocity 
tends to a constant value as C increases. 
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FIG. 5. Mean rate of energy change versus velocity, for 
3D simulations with grid spacing 0.234 (256 x 64 x 64 points). 
Parameters are C = 1000, e = 0.0625, rj = 1, a = 3.0, Wo = 1, 
and A ~ 0.281. The speed of sound at the condensate center 
c s ~ 2.54 is represented by the dashed line. For U b = 40 
(bullets) the critical velocity is v c ~ 0.13c s , while for U a h = 20 
(triangles) it is v c ~ 0.20c s . 

Simulations in 3D were performed, and the mean en- 
ergy transfer rate as a function of velocity is presented in 
Fig. H Similar behaviour to 2D is observed, with smaller 
critical velocities: a result of the beam intersecting the 
condensate edge where the speed of sound is lower. Ac- 
cordingly, vortex lines first appear in these regions and 
penetrate into the center. This conclusion agrees with 
the experiment ]7|Jlq|, where a relatively low critical ve- 
locity (v c ~ 0.26c s jwas measured. The dependence of v c 
on U a h and C was found to be similar to 2D, where e.g. 
v c ~ 0.29c s for C = 4000 and U b = 35. Enhanced heat- 
ing is also observed for v c > c s , due to phonon emission 
between the extrema. 

In this paper, we have studied the role of vortex for- 
mation in the breakdown of superfluidity, by an oscillat- 
ing object in a trapped Bose-Einstein condensate. We 
find that at low object velocities, energy is transferred 
by phonon emission at the motion extrema, while a much 
larger energy is transferred above the critical velocity for 
vortex formation. To generalize these conclusions to re- 
alistic experimental situations, the model should include 
the non-condensed thermal cloud. Energy would then be 
transferred from the condensed to the thermal cloud by 
phonon damping Q, or vortex decay p^j] . 
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